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A combined kinetic and computational study on our tryptophan-based bifunctional thiourea catalyzed 
asymmetric Mannich reactions reveals an apparent negative activation enthalpy. The formation of the 
pre-transition state complex has been unambiguously confirmed and these observations provide an 
experimental support for the formation of multiple hydrogen bonding network between the substrates and 
the catalyst. Such interactions allow the creation of a binding cavity, a key factor to install high 
enantioselectivity. 

Hydrogen bonding is an essential interaction in biological systems 1 . In addition to their role in structural 
determination, it is crucial in electrophile activation in biocatalysis 2 ' 3 . For instance, amide hydrolysis 
catalyzed by serine protease provides a useful starting point for the discussion of electrophile activation 
by hydrogen bonding interaction in enzymatic catalysis 4 . The mechanistic study of [3,3]-sigmatropic rearrange- 
ment catalyzed by chorismate mutase has also demonstrated hydrogen bonding interaction's importance in 
binding the substrates in a preferred conformation 5 " 8 . Research in organocatalysis has sought to develop small 
molecules to utilize hydrogen bonding interactions to achieve enzyme-like activity and selectivity 9 " 14 . While 
numerous successful examples have been demonstrated in the past two decades, kinetic understanding and 
evidence of these systems is still limited. 

Thiourea derivatives have been regarded as effective Lewis acid-like catalysts for their ability in electrophilic 
activation through hydrogen bonding 9 . The enantioselective hydrocyanation of imines (Strecker reaction) cat- 
alyzed by Schiff base derived chiral thiourea catalysts 15 and subsequent kinetic and computational studies 16 by 
Jacobsen and co-workers marked one of the early examples of utilizing this form of activation in asymmetric 
catalysis. The first example of bifunctional thiourea catalyst bearing tertiary amine was reported by Takemoto and 
co-workers, and they showed that such catalysts led to excellent enantionselectivity in the Michael reaction of 
malonate esters with nitroalkenes 1718 . The activation model proposed indicated that the nitroolefin interacts with 
the thiourea moiety via double hydrogen bonding, while the tertiary amine holds the malonate via single 
hydrogen bonding in the activated enol form. In the past decade, novel chiral thiourea derivatives have been 
developed and proliferated as the most prominent hydrogen bond donor catalysts in a wide variety of asymmetric 
organic reactions 10 ' 16 ' 19 " 23 . 

We recently demonstrated that tryptophan-based bifunctional thiourea derivatives are capable of efficiently 
promoting the asymmetric Mannich reaction of fluorinated ketoesters to afford highly optically enriched flu- 
orine-containing molecules with adjacent quaternary and tertiary stereocenters (Figure l) 24 . High enantioselec- 
tivity and diastereoselectivity were observed with a wide range of aromatic and heteroaromatic ketoesters with an 
optimal bifunctional thiourea-based organocatalyst. Preliminary computational studies also suggested the 
involvement of a "pre-transition- state complex" with multiple hydrogen bonding interactions, and in turn could 
be responsible for the excellent stereochemical control. In a bid to further deepen our mechanistic understanding 
of such a catalytic system, we attempted kinetic studies in which experimental data revealed an uncommon 
apparent negative activation enthalpy. Even though there is precedence of experimentally observed negative 
activation enthalpy in reactions 25 " 29 , the phenomenon in organocatalytic systems is nascent revealing insightful 
chemical behavior. Concomitant computational studies supported experimental findings and from a theoretical 
viewpoint we were able to elucidate the origin of the stereochemical outcome and to deduce that the observed 
apparent negative activation enthalpy arose from a highly extensive and favorable H -bonding effect. 
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Figure 1 | Mannich reaction of fluorinated ketoester and imine. 



Results 

Determining the overall rate of reaction. Kinetic studies were 
conducted to verify the termolecularity of the reaction (Supporting 
information). When the reaction was carried out with 10 equivalent 
excess of ketoester B, the -ln([A]/[A] 0 ) plot versus time of imine A 
gave a straight line (R 2 = 0.999), which indicated that the reaction is 
first-order with respect to A. In the same light, it was determined that 
the reaction is first- order with respect to B with 10 equivalent excess 
of imine A. The order in catalyst was also examined by plotting the 
kinetic rate constant (fc 0 & s ) against the loading of the catalyst C (Trp- 
1), which satisfy that the reaction is first-order with respect to C. As 
the overall reaction rate order is 3, this confirms that this catalyst 
system is a termolecular process and reaction rate would be 
dependent on the concentrations of A, B and C. Another general 
base catalyst, 1,4-diazabicyclo [2.2.2] octane (DABCO) was used for 
the comparison with C. Employing previous protocol, it was found 
that the reaction rate still first-order with respect to each substrate 
and the catalyst. 

Apparent negative activation enthalpy. Activation parameters were 
acquired by the means of Eyring plot analysis with the measurement 
of reactions rates at 10-40°C 30 . In sharp contrast to analogous studies 
done by Takemoto 17 , it was found that the reaction rate decreased as 
the reaction temperature was raised within the experimentally 
accessible temperature range (Table 1) and the activation enthalpy 
AH obs x was extrapolated to be negative (Figure 2)! The activation 
entropy AS ofe J of —74.7 e.u. agrees well with a termolecular reac- 
tion. Large activation entropy observed is not without precedence as 
reported in other pioneering computational and kinetic works in 
thiourea catalysis 26,27 . It was noted that the corresponding DABCO 
studies found the observed activation enthalpy to be positive and the 
reaction rates increased with elevated temperatures (see supporting 
information). 



Table 1 | Activation parameters for Mannich reaction of imine 
and fluorinated ketoester catalyzed by thiourea C at different 
temperatures 
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DFT investigations on Mannich reactions with Trp-1. This 
reaction was modeled with DFT calculations to elucidate the 
observed enantioselectivity and to appraise the energetics. It is 
well -accepted that the pendent amine on the catalyst can first act 
as a base to assist the formation of the ketoenolate from the P- 
ketoester substrate. The resulting protonated amino group together 
with the thiourea moiety can then bind the ketoenolate and activate 
the imine electrophile. Two pre-transition-state models of different 
reactant arrangements derived from theoretical investigations have 
been proposed in the literature (Figure 3) 31 " 40 . Our initial efforts were 
thus focused on the identification of the pre-transition-state complex 
between the substrate and catalyst. Hydrogen bonding is believed to 
be the keystone interaction for the binding of catalyst and substrates 
leading to a spatially preferred conformation dictated by the catalyst, 
a phenomenon that is a reminiscence of enzyme-substrate(s) 
complexation in biological systems. 

Results from our modeling indicate that the computed pre-trans- 
ition-state structure pTS-I (model A; +19.4 kcal mol" 1 relative to 
starting materials) is similar in energy to pTS-II (model B; 
+ 19.0 kcal mol -1 ), but the transition state energy barrier of TS-I 
via pTS-I was found to be significantly lower than that of TS-II via 
pTS-II by 14.0 kcal mol" 1 . The vast difference in activation barrier 
energies for TS-I and TS-II could be rationalized that both electro- 
phile and nucleophile are bound closer and better spatially fitted in 
the catalyst's 'assembly pocket' in pTS-I as compared to pTS-II (C-C 
bond distance for pTS-I is 3.906 A and pTS-II is 5.303 A). In other 
words, a more transition state-like orientation of the substrates was 
realized in pTS-I. Moreover, for pTS-I, the oc-C-H bonds (to the 
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Figure 2 | Eyring plots showing temperature dependence behaviour of 
the Mannich reaction catalyzed by Trp-1 measured by 1 !! NMR. 
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Figure 3 | Proposed models for catalyst- substrate complex. 



ammonium group and to the indole ring) are in proximity (2.1 to 2.4 
A) with the ester carbonyl oxygen (Figure 4), alluding to the pos- 
sibility of the formation of a non-classical C-H s+ ...O s ~ binding 
pocket. Such secondary interactions have been observed to play an 
important role in molecular recognition and stereoselective catalytic 
processes 41 " 43 . To confirm that such short bond distances are result of 
hydrogen bonding interactions rather than simple steric arrange- 
ment, two-center Mayer bond order analysis 44 was carried out. The 
bond indices were calculated from canonical molecular orbitals in 
the atomic orbital basis and can be used to analyze the covalent 
bonding between atoms. The computed bond orders of 0.057 and 
0.032 for these a-C-H O bonding interactions are more significant 
with reference from our recent DFT studies with different derivatives 
of thiourea for catalytic organic transformation 45,46 . The N-H O 



TS-I 

+29.1 kcal/mol 
-1.1 kcal/mol 




TS IV 

+30.4 kcal/mol 
+0.3 kcal/mol 



Figure 4 | Geometry of key intermediate and transition states showing bond distances (A) of non-covalent interactions (blue) (mPWlPW9 1/6-31 + 
G(d,p)). Mayer bond order for pTS-I are in parenthesis and blue. Reported energies are relative to starting materials: free energy in black non-italic 
bold and enthalpy in black italic bold. Dotted lines show hydrogen bonding (black) and C-C bond formation (orange). Non-interacting hydrogen atoms 
have been omitted for visual clarity. Atoms are color coded for each element: carbon (grey); hydrogen (white); oxygen (red); nitrogen (blue); fluorine 
(green); sulfur (yellow). Insert showing the Newman projection of the transition states from the perspective of the newly formed C-C bonds. 
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hydrogen bonding interactions showed similar to prominently 
higher bond orders of 0.040, 0.098, and 0.155. Both the C-H and 
N-H hydrogen bonding interactions contribute to the creation of the 
binding cavity and subsequently prompt the ketoenolate to be bound 
such that providing overall stabilization of the charged ketoenolate 
and locking the nucleophile more rigidly. The imine on the other 
hand binds via hydrogen bonding to the pendent ammonium group, 
which firstly enhances the electrophilicity of the sp 2 carbon and 
secondly, fits well above the ketoenolate and is predisposed to 
nucleophilic addition for the C-C bond formation (Figure 4, top). 

The imine molecule could bind in two different faces relative to the 
ketoenolate: Re or Si (as viewed by the approach of the nucleophile) 
and the ketoenolate could also bind to the protonated catalyst in two 
different orientations, thus determining the overall stereoselective 
outcome. Computationally, it was found that the transition state 
leading to the .Re-attack, TS-I, is 15.0 kcal mol" 1 more favorable than 
the S/-attack through TS-III. The preference for TS-I is obvious as 
the s-trans-imine can be bound more closely to the ketoenolate. In 
order for the S/-face attack on the imine to materialize, the imine has 
to be the s-cis form which is energetically penalized due to the steric 
repulsion between the Boc and the 4-methylphenyl groups (TS-III, 
Figure 4). Modeling for the second orientation of the ketoenolate 
binding to the catalyst was carried out and the corresponding trans- 
ition state TS-IV is slightly endergonic by 1.3 kcal mol" 1 compared 
to TS-I, accounting for the 9 : 1 diastereoselectivity for the reaction. 
These observations provide the rationale to the constant excellent 
enantioselectivity obtained for all ketoester substrates yet moderate 
diastereoselectivity when the ketoester has substituents of similar size 
(e.g. methyl and ethyl groups) 24 . It is also noted that the calculated 
differences in activation free energies of these transition states were 
similar to those in activation enthalpies, suggesting that the observed 
selectivities are not due to the entropic control in Sohtome and 
Nagasawa's systems 42,47 . 

When measuring the reaction rate, one is looking at the overall free 
energy barrier, i.e. TS-I, and it occurs that the overall enthalpy 
change of transition state TS-I is — 1.1 kcal mol" 1 relative to starting 



substrates. Comparing the enthalpy change for intermediate pTS-I to 
starting materials ( — 7.2 kcal mol" 1 ) suggests a distinct feature of an 
exothermic elementary step due to formation of hydrogen bonds 
with C and substrates. As our experiments revealed inverse depend- 
ence of reaction rate with temperature, it is clear from the reaction 
profile how negative activation enthalpy originated (Figure 5). 

Discussion 

From our combined theoretical and kinetic mechanistic study for the 
enantioselective bifunctional thiourea catalyzed Mannich reaction, 
first experimental evidence was revealed for an apparent negative 
activation enthalpy. Since the apparent enthalpy of activation is the 
sum of enthalpies of the complexation in pre-transition state and of 
the bond formation step, this important observation suggests that the 
multitude of non-covalent bonding interactions between the catalyst 
and substrates certainly translate to a more exothermic enthalpic 
change for the complexation step (pTS-1, Figure 5) with AH = 
— 7.2 kcal mol" 1 . As a result of the highly exothermic adduct-form- 
ing step and only slightly endothermic ( + 6.1 kcal.mol) correspond- 
ing activation enthalpy for the elementary C-C bond formation step, 
the observed enthalpy was analyzed to be —3.1 kcal mol" 1 (AH^*). 
The empirical treatment of AG 0 b s t = AH obs * — TAS^* for our 
kinetic measurements, for which /SH ohs l and AS 0 & S * are assumed to 
be constant for a small change in temperature, allowed us to observe 
an apparent negative activation enthalpy due to the exothermic 
change from starting intermediates to the transition state as being 
clarified computationally at the same time. These observations un- 
ambiguously confirm the multistep nature of the mechanism in the 
chiral bifunctional thiourea catalyzed reactions, i.e. a pre-transition- 
state complex is formed. Consistent with the kinetic indications, our 
calculation results further suggest the formation of multiple hydro- 
gen bonding interactions between the substrates and the catalyst 
after the pendent amine group deprotonates the pro-nucleophile 
(ketoester), supporting the analogous model proposed by Papai 
and co-workers 48 . These interactions allow the creation of a binding 
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cavity to direct the incoming imine substrate to install high enantios- 
electivities. 

Methods 

Kinetics experiments were recorded on a 500 MHz spectrometer. Chemical shifts 
were reported in parts per million (ppm), and the residual proton peak was used as an 
internal reference: proton (in Benzene-d 6 : § 7.16). Benzene-d 6 was dried from 
sodium/benzophenone. Stock solutions in deuterated benzene were used in kinetics 
experiments and 1,4-dimethoxybenzene was used as internal standard. 
Computations were carried out with density functional theory (DFT) employing 
Gaussian 09 suite of program 49 . The modified Perdew-Wang exchange functional 
mPWlPW91 theory 50-52 and the all electron split- valence Pople basis set 6— 31 + 
G(d,p) 53-55 containing polarization functions on heavy atoms and hydrogen and 
diffuse functions for heavy atoms was used as it was shown to be accurate for mod- 
eling Mannich reactions 56 . Frequency analyses were carried out to verify minimum 
structures having only positive eigenvalues of the Hessian matrix or transition state 
structures showing a single negative eigenvalue. 
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